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Kinetic  Theory  of  Field  Evaporation  of  Metals 
L.C.  Wang  and  H.J.  Kreuzer 
Department  of  Physics,  Dalhousie  University 
Halifax,  N.S.  B3H  3J5,  Canada 


Abstract  -  Starting  from  adiabatic  energy  curves,  calculated  as  a  function  of  elec¬ 
tric  field  strength  we  develop  a  perturbational  method  to  construct  diabatic  states 
which  form  the  basis  to  compute  the  temperature  dependent  ionization  probabilities 
for  field  evaporation  of  metals.  Employing  a  master  equation  we  calculate  the 
energy  dependent  ion  yield  as  a  function  of  field  strength  and  temperature,  and 
extract  the  field  dependence  of  the  activation  barrier  and  the  prefactor.  As  a  func¬ 
tion  of  field  strength  the  activation  barrier  for  field  evaporation  decreases  mono- 
tonically  for  tungsten,  whereas  the  prefactor  increases  up  to  a  field  strength  of 
4.5V/S  for  tungsten  due  to  an  increase  in  the  ionization  probability  and  then  decre¬ 
ases  due  to  changes  in  the  surface  potential. 
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1 .  Introduction 

Field  evaporation  is  the  removal  of  lattice  atoms  from  the  surface  of  a  field 
emitter  tip  as  singly  or  multiply  charged  positive  ions  in  an  electric  field  of  the 
order  of  volts  per  angstroms  [1,2].  The  term  field  desorption  is  usually  reserved 
for  the  process  of  removing  field-adsorbed  atoms  or  molecules  from  field  emitter 
tips  [3-5].  Field  evaporation  and  field  desorption,  crucial  in  cleaning  and  preparing 
field  emitter  tips,  are  thermally  activated  processes;  as  such  their  rate  constants 
can  be  parametrized  as 

rd  =  av  exp(-Q(F)/kBT)  (1) 

Here  Q(F)  is  the  field  dependent  height  of  the  activation  barrier  to  be  overcome  by 
the  desorbing  particle.  We  have  split  the  prefactor  into  an  attempt  frequency  v  and 
an  accommodation  coefficient  a,  both  being  field  erd  also  weakly  temperature  depen¬ 
dent.  The  minimum  field  strength  beyond  which  the  activation  barrier  vanishes,  thus 
leading  to  field  evapration  at  low  temperatures,  is  called  the  evaporation  field 
strength;  it  varies  from  2.5V/H  for  Ti  to  6.1V/JI  for  W  with  a  typical  experimental 
error  margin  of  10-20  $  [5].  Ernst  has  measured  Q(F)  and  av(F)  for  Rh  [6], and  Kel¬ 
logg  [7]  has  presented  data  for  W  in  the  field  range  4.7-5.9V/iL 

Two  phenomenological  models  have  been  proposed  to  calculate  the  activation  energy 
Q(F) .  In  the  "image-force"  model  [2]  field  evaporation  is  envisaged  as  the  activation 
of  an  ion  of  charge  ne  over  an  activation  barrier  that  results  from  the  superposition 
of  the  field  potential  -neFz  (assuming  a  constant  electric  field),  and  the  image  pot¬ 
ential  of  the  ion  -neAz.  In  the  "charge-exchange"  model  [8]  one  assumes  that  ioni¬ 
zation  and  desorption  occurs  at  the  crossover  point  between  the  atomic  and  ionic 
(diabatic)  potential  energy  curves.  A  confrontation  of  these  models  with  experimen¬ 
tal  data  has  been  presented  by  Kellogg  [7];  Forbes  [9]  has  presented  several  criti¬ 
cal  assessments.  An  early  microscopic  calculation  of  Q(F)  by  Kahn  and  Ying  [10]  was 
based  on  the  local  density  approximation  of  the  density  functional  theory,  treating 
the  metal  as  a  jellium.  Kingham  [11]  has  presented  some  preliminary  results  for  the 
field  evapration  of  W  obtained  within  a  tight  binding  cluster  model.  Kreuzer  and 
Nath  [12]  have  presented  a  microscopic  theory  of  field  evaporation  in  which  the 
electronic  properties  of  the  metal  are  calculated  within  a  tight  binding  cluster 
approach  based  on  the  ASED-MO  method  [13]  with  local  electric  fields  taken  from 
selfconsistent  jellium  calculations  [1 4,15].  They  extract  the  activation  energy  Q(F) 
from  adiabatic  ground  state  energy  curves,  and  find  a  scaling  law  which  predicts 
evaporation  field  strengths  within  10$ -20?  of  experimental  values. 

To  extend  such  ground  state  energy  calculations  to  a  kinetic  theory  of  field 
desorption  and  evaporation,  Kreuzer  et  al.  [16-18]  have  employed  a  master  equation 
to  calculate  the  energy  dependent  ion  yield  in  thermal  field  desorption  as  a  function 
of  field  strength  and  temperature,  applying  their  method  to  thermal  field  desorption 
of  helium.  In  this  paper  we  will  briefly  review  this  theory  and  then  use  it  to  study 
field  evaporation  of  metals. 

In  the  next  section  we  will  review  the  calculation  of  adiabatic  and  diabatic 
energy  curves  in  the  presence  of  high  electriu  fields  at  surfaces.  The  latter  will 
then  serve  as  the  framework  to  compute  the  temperature  dependent  ionization  and 
neutralization  probabilities  of  atoms  on  the  surface  of  a  metal.  In  the  last  section 
we  will  present  numerical  data  for  field  evaporation  of  W.  To  summarize  our  res¬ 
ults:  we  find,  as  in  our  earlier  study  [12],  that  the  activation  barrier  against  field 
evaporation  decreases  (monotonically  for  most  metals,  such  as  tungsten,  for  which  we 
find  good  agreement  with  Kellogg's  data  [7])  whereas  the  prefactor  initially  incre¬ 
ases  (up  to  i).5V/A  for  tungsten)  due  to  an  increase  in  the  ionization  probability,  and 
then  decreases  due  to  changes  in  the  surface  potential  of  the  desorbing  atom. 
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2.  Adiabatic  and  diabatic  states 

Field  evaporation  takes  place  predominantly  at  kink  and  step  sites  due  to  the  fact 
that  electric  fields  are  enhanced  at  such  sites.  To  set  up  our  theory  we  rather 
consider  a  situation  where  a  single  metal  atom  is  present  on  top  an  otherwise  per¬ 
fect  crystal  plane,  a  distance  R  away  from  the  topmost  ion  cores.  We  will  refer  to 
this  isolated  atom  as  the  adatom.  The  hamiltonian  of  the  system  can  be  written  as 

H  =  Tn  +  He(rltr2,...;R)  (2) 


where 


TN  =  32/9R2 

is  the  kinetic  energy  of  the  atomic  nucleus  and 


(3) 
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is  the  hamiltonian  of  the  electrons  at  positions  r1,r2,...  Ve  includes  the  Coulomb 
interactions  between  the  electrons,  between  the  electrons  and  the  nuclei,  and 
between  all  nuclei  (metallic  and  atomic). 

Let  us  fix  the  position  of  the  adatom,  thus  setting  its  kinetic  energy  (3)  equal 
to  zero.  Physically  this  implies  that  the  electronic  degrees  of  freedom  follow  the 
nuclear  motion  instantly.  We  can  then  diagonalize  He  (in  practice,  after  approximat¬ 
ing  it,  e.g.,  by  a  tight  binding  hamiltonian  or  using  density  functional  theory)  to 
obtain 


He(r  ,,r2, . .  .;R)  £^(r  i,r2, .  •  .;R)  —  V^(R)  Ci(r  j,r2,...;R)  (5) 

where  the  are  adiabatic  many  electron  wavefunctions.  The  lowest  eigenvalue  of 
(5),  V0(R)  represents  the  ground  state  of  the  system  and  corresponds  to  the  adiabatic 
binding  energy  curve.  Lifting  an  electron  from  the  highest  occupied  level  (in  the 
ground  state)  to  the  lowest  unoccupied  one,  generates  the  energy  curve  (or  rather 
surface,  because  R  is  a  three-dimensional  space  coordinate)  of  the  first  excited 
state  etc. 

To  take  the  nuclear  motion  of  the  adatom  into  account  we  now  proceed  to  con¬ 
struct  diabatic  states.  Returning  to  the  hamiltonian  (2)  we  try  to  diagonalize  it  by 
solving  Schrodinger's  equation 

H  ¥<*  ‘  Ea  Ya  (6) 

via  an  expansion 


l*,a^r i»r2,.«.;R)  =  Y  Ci(ri,r2,...;R)  xia(R)-  (?) 

i 

Inserting  (7)  into  (6),  multiplying  with  Ci*  and  integrating  over  the  electronic  deg¬ 
rees  of  freedom,  r,,  r2,....,  we  obtain  a  coupled  set  of  equations 
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are  the  matrix  elements  of  the  first  and  second  order  differential  operators.  Put¬ 
ting  the  right  hand  side  equal  to  zero  results  in  the  Born-Oppenheimer  approximation 
according  to  which  the  nucleus  just  moves  in  the  i-th  adiabatic  potential. 

To  solve  (8)  amounts  to  obtaining  an  exact  solution  to  the  problem  by  diagonalizing 
an  infinite  matrix.  To  reduce  the  task  to  manageable  proportions  one  hopes  that  in  a 
given  problem  only  a  limited  number  of  adiabatic  states  contribute  to  the  set  of 
diabatic  states  involved  in  the  process  to  be  studied.  E.g.  in  field  desorption  and 
field  evaporation,  only  the  lowest  diabatic  states  for  a  neutral  atom  and  for  ions 
seem  important.  They  can  thus  be  obtained  by  a  unitary  transformation  from  the  cor¬ 
responding  adiabatic  states.  The  latter  being  known  explicitly  from  the  tight  bind¬ 
ing  or  density  functional  calculations  alluded  to  above,  the  problem  is  solved  in 
practice. 

Let  us  introduce  new  states 

x  =  A  n  (ID 

where  xT=(XnXz . )  is  the  transpose  of  the  column  vector  x  containing  the  adiabatic 

states  as  its  components  and  A  is  a  unitary  matrix  which  is  an  explicit  function  of 
R.  To  obtain  a  diabatic  solution  of  (8)  we  choose  the  transformation  A  such  that  the 
first  order  derivative  term  vanishes,  i.e.  we  impose  the  condition  [19] 

(—  +  T^1))  A  =  0  (12) 

We  thus  get  from  (8) 

(-  fft  -  W  +  -(R)  "  "  =  0  ( 1  3) 

where  n  are  diabatic  states.  In  (13)  I  is  a  unit  matrix,  E  is  diagonal  and  the  dia¬ 
batic  interaction  matrix,  no  longer  diagonal,  is  given  by 

W  =  AfV  A  (1 10 

in  terms  of  the  adiabatic  energy  curves  defined  in  (5).  Dropping  the  coupling  terms 
Wjj  from  (13)  determines  uncoupled  diabatic  states  subject  to 

(-  Jp  +  Wii(R)  -  Eiv)  n°iv  = 


0 


(15) 
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Thus  n°iv  is  the  vth  eigenstate  with  energy  Eiv  of  the  ith  diabatic  potential  W^R). 
E.g.  W00  would  be  the  potential  for  a  neutral  atom  approaching  a  surface,  and  Wn 
that  of  an  ion.  The  off-diagonal  terms  Wjj  couple  these  states  together;  in  the 
above  example,  W01  is  responsible  for  ionization  or  neutralization. 

Combining  (7)  and  (11),  we  can  rewrite  the  total  wavef unction  as 

T(r;R)  =  eT(r;R)X(R) 

=  ?T(r;R)A(R)n(R) 

=  £T(r;R)n(R)  (16) 

introducing  the  diabatic  many  electron  wavefunctions,  £T(r;R)=( £0,  .. .) . 

Before  we  continue  with  the  explicit  calculation  of  adiabatic  and  diabatic  states, 
we  want  to  briefly  discuss  the  meaning  of  these  states.  To  this  end  we  assume  that 
in  the  process  of  field  desorption  or  field  evaporation  only  two  levels  participate 
with  the  lower  electronic  states  remaining  unchanged  as  spectators  throughout  the 
process.  We  denote  by  $a  the  highest  occupied  single  electron  level  with  energy  ea 
in  the  isolated  adatom  and  by  the  lowest  unoccupied  level  of  the  same  symmetry 
with  energy  em  above  the  Fermi  energy  of  the  isolated  metal.  In  the  interacting 

system  of  adatom  plus  metal  these  two  levels  combine  to  two  molecular  orbitals  <t>g 

and  4>e  of  lower  and  higher  energy,  respectively,  see  Fig.1.  Assuming  that  we  can 
approximate  the  adiabatic  many  electron  wavefunctions  as  Slater  determinants,  we 
find  that  the  adiabatic  ground  state  wavefunction  for  N  electrons 

C0(r;R)  =  |$14>2...<J>N-i<5>g|  (17) 

is  composed  of  the  spectator  states,  4>i  for  i=1,...N-1,  and  <t>g,  whereas  the  lowest 
excited  many  electron  wavefunction 

Ci(rsR)  =  |M2...<tN_i<s>e|  H8) 

has  <t>g  replaced  by  <t>e.  Assuming  that  with  the  adatom  close  to  the  metal  ea  is  lower 

than  em,  <t>g  and  4>e  are  dominated  by  <f>a  and  $m,  respectively.  Thus  £0  represents  a 

many  electron  state  in  which  both  adatom  and  metal  are  in  their  respective  neutral 
ground  states,  whereas  s,  describes  a  situation  with  an  electron  removed  from  the 
adatom  and  put  onto  the  metal. 

We  next  look  at  a  situation  where,  in  the  presence  of  an  electric  field  the 

adatom  is  so  far  away  from  the  surface  that  ea  is  higher  than  em.  As  a  result  the 

molecular  orbital  <J>g  is  now  mainly  of  4>m  character,  whereas  4>e  is  mostly  <j>a,  imply¬ 
ing  that  £0  now  describes  a  system  with  the  adatom  being  ionic  and  the  removed  elec¬ 
tron  in  the  metal  above  the  Fermi  energy,  whereas  £i  represents  neutral  adatom  and 
metal . 

In  our  previous  paper  [18]  we  had  shown  that  for  two  relevant  states  the  trans¬ 
formation  A_  in  (11)  can  be  represented  by  a  unitary  2x2  matrix  which  can  be  written, 
quite  generally,  as 

.  [  cosQ(R)  sinO(R)"]  . 

L  -sme(R)  cosQ(R) J  •  uy; 

Noting  that  the  matrix  t(^  is  antisymmetric,  we  reduce  the  matrix  equation  (12)  to  a 
simple  vector  equation 
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which  can  be  solved  as  a  line  integral 


rR 


OCR)  =  - 


dR' •<  £0 


(21) 


Having  0(Ro)=O  implies  that  A(R0)=I_  and  thus  that  at  R0  diabatic  and  adiabatic  states 
coincide  which  is  the  case  far  from  the  surface.  The  diabatic  interaction  matrix 
(1*0  is  now  given  explicitly  by  its  components 

W00(R)  =  cos20(R)  V„(R)  +  sin2Q(R)  VX(R)  (22) 

Wn(R)  =  cos20(R)  V,(R)  +  sin20(R)  V0(R)  (23) 

Wol(R)  =  4  sin20(R)  [V0(R)  -  VL(R)]  .  (2h) 


To  evaluate  (21),  we  observed  that  in  field  desorption  adatoms  will  leave  the  sur¬ 
face  along  the  steepest  field  gradient,  i.e.  perpendicular  to  the  surface  so  that  we 
can  neglect  any  lateral  interactions  and  O=0(z)  depends  on  the  distance  z  from  the 
metal  only  so  that 


0(z) 


z 

dz/dr  <tg(r;z)  3/3z  <fe(r;z) 

00 


(25) 


provided  that  the  adiabatic  many  electron  wavefunctions  are  given  by  Slater  determi¬ 
nants.  This  completes  the  construction  of  diabatic  states  for  situations  where  two 
states  are  important.  We  should  note,  however,  that  even  in  the  case  of  thermal 
field  desorption  of  helium,  this  is  not  strictly  true.  Recall  that  field  ionization 
of  the  adatom  can  take  place  when  its  highest  occupied  level  rises  above  the  lowest 
unoccupied  level  in  the  metal.  For  a  given  field  strength,  F,  this  happens,  if  the 
adtom  is  at  a  distance  zc(F).  If  the  desorbing  atom  is  at  distances  z  larger  than 
zc,  the  electron  can  tunnel  into  higher  unoccupied  levels  of  the  metal.  Thus  for 
given  F  and  z,  only  two  states  participate  in  the  ionization  process  and  thus  in  the 
construction  of  diabatic  states,  although  for  the  overall  field  desorption  process  a 
band  of  metal  states  is  relevant.  Note,  however  that,  because  the  tunneling  proba¬ 
bility  decreases  rapidly  as  a  function  of  distance,  this  band  of  tunneling  states  is 
rather  narrow.  Mimicking  the  metal  by  a  finite  cluster  of  metal  atoms  obviously 
does  not  produce  a  band  structure  but  only  a  set  of  discrete  levels  which  one  has  to 
broaden  with  a  width  corresponding  to  the  width  of  the  ionization  zone 

Dealing  with  thermal  field  evaporation  of  metal  ions,  the  situation  becomes  more 
complex  because  it  is  less  evident  how  to  assign  molecular  orbitals  of  the  metal 
plus  adatom  complex  to  either  the  adatom  or  the  metal,  particularly  if  the  latter  is 
approximated  by  a  finite  cluster  itself.  Many  states  then  become  important  corres¬ 
ponding  to  different  stages  of  ionization  and,  more  importantly,  of  different  outgo¬ 
ing  momentum.  In  principle,  we  can  use  the  above  formalism  resulting  in  the  con¬ 
struction  of  large  matrices  A.  Instead,  we  will  advance  a  perturbative  approach  to 
the  construction  of  diabatic  states  that  starts  from  the  observation  that  the 


7 


transfer  of  an  electron  from  the  adatom  to  the  metal  takes  place  when  the  highest 
occupied  level  on  the  adatom  rises  above  the  Fermi  level  of  the  metal.  If  this  hap¬ 
pens  with  the  adatom  at  position  zc,  we  can  linearize  the  energy  of  the  highest 
occupied  level  as 

ea(z)  =  em  +  Ae(z)  (26) 

with  Ae(zc)  =0.  We  further  write  for  the  molecular  orbitals,  restricting  ourselves 
to  two  levels  for  the  moment, 


^g  "  cgm<t>m  +  cga^a 
=  cem<t’m  +  cea<Jla 
and  get  to  second  order 


_L 

i-(f-g)s 

✓2 

/l+f2-fg 

1 

f-g 

/2 

J 1  +f  4-fg 

1 

1-(f+g)S 

/2 

/ 1  +f  2+g 

c  -  JL  f  +  g 

gm  “  . - — 

*  2  J 1  +f  2+g 


where 


f(z)  =  Ae(z)/(2(V-emS) 


(27) 


(28) 


g  =  /T+F5 


V  -  <  1  h  |  <t>a  > 

S  =  < 0m  I  4>a >  (29) 

The  transformation  angle  (25)  is  then  given  by 

0(z)  =  4tg"*f(z)  -  it/U  (30) 

where  h  is  the  single  electron  hamiltonian.  Having  0(z),  we  can  next  calculate  the 
diabatic  potentials  Wjj  from  (22-24).  In  particular,  note  that 

Vo-V»  =  -g(V-emS)  (3D 


so  that 
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W01(z)  =  V(z)  -  emS(z)  (32) 

The  diabatic  many  electron  wave  functions  are  then  given  by 

£0(r;R)  =  (cgmcos0  -  cemsin0)  |4>,4>2...<tN_1 4>m  | 

+  (cgacos0  -  ceasin0)  4>a  | 

5j(r;R)  =  (cgmsin0  +  cemcos0)  | <*>14>2 — «*>N-1  4>m  | 

+  (cgasin0  +  ceacos0)  |4>,4>2...$n_i  4>a  |  (33) 

where,  with  (30), 

sin0(z)  =  -/I -f/g 

cos0(z)  =  /T+f7g  (3*0 

so  that  (33)  reduces,  properly  normalized,  to 

Ur;R)  =  (1-S2rl/2  (|*1*a...*N_1*a|  "s  <!>m  |) 

6i(r;R)  =  |Ma**.*N-l*m|  (35) 

where  S  is  the  overlap  of  levels  <f>a  and  <j>m,  as  given  in  (29). 

To  assess  this  perturbative  approach  we  have  recalculated  a  few  numbers  relevant 
for  field  desorption  of  helium  from  tungsten.  E.g.,  in  a  field  of  8.5V/i( 
woi(zc)=1 «96meV  using  the  exact  transformation  and  1.78meV  employing  the  approximate 
formula  (32).  In  the  relevant  interval  of  z  around  zc  the  error  is  never  larger  than 
30%. 

We  now  want  to  generalize  the  perturbative  approach  to  situations  where  several, 
or  many  states  on  the  adatom  and  in  the  metal  participate  in  the  ionization  process. 
We  first  rewrite  the  interaction  matrix  of  the  hamiltonian  in  the  diabatic  basis  as 

W(R )  =  /dr  _fj(r;R)  He(r;R)  ST(r;R)  (35) 

and  ass  :me  again  that  the  many  electron  wave  functions  can  be  approximated  by 
single  Slater  determinants.  We  denote  by  <pa  the  wave  function  of  the  highest  occu 
pied  level  on  the  isolated  atom  and  by  4>m  ^  with  k=l,2,...  single  electron  wave  func¬ 
tions  of  energy  em  k  above  the  Fermi  level  of  the  isolated  metal.  To  ensure  ortho¬ 
gonality,  we  define 

|$a>  =  NaPm|*a>  (36> 

where  the  projection  operator  is  given  by 


1 4>m ,k><^m ,k 
k 


(37) 


so  that  the  normalization  constant  becomes 


-2 


(38) 
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N 


a 


<  4>m  ,k  I  l^a>' 
k 


The  diabatic  many  electron  wave  function  of  lowest  energy  with  both  adatom  and 
metal  neutral  is  then  given  by 

So(r;R)  =  |*1*a...*N_i*a|  (39) 

The  excited  states  are  described  by 

Ck(r;R)  =  |<M2...<t'N-1<t>mfk|  (^0) 

having  one  electron  transferred  from  the  adatom  t  >  the  k-t*  level  above  the  Fermi 
energy  of  the  metal.  Writing  the  electronic  hamiltonian  as  a  sum  of  single  electron 
contributions 


N 

H e Cr* ; R )  =  h(ri;R) 
i=1 

and  inserting  (39)  and  (40)  into  (35)  we  get  for  the  offdiagonal  terms 
=  Na(Vk  -  emikSk) 


(41) 


(41) 


where 

^k  =  < ^a I ^  I'J’m  .k^  (42) 

is  the  interaction  between  states  $a  and  4>m  k,  and 

^k  =  < 't’a  I  ,k>  (43) 

is  the  overlap  between  them.  This  completes  the  perturbative  approach  tc  the  calcu¬ 
lation  of  the  diabatic  interaction  matrix  (35). 

5.  Results 

In  our  previous  paper  [18]  we  have  outlined  the  derivation  of  kinetic  equations  that 
control  field  ionization,  field  desorption  and  field  evaporation.  In  particular,  we 
derived  an  expression  for  the  yield  of  singly  charged  ions 


Yion  =  ]Lt0+(u,v)  exp(-E0M/kBT)/ ^  exp(-E0u/kBT)  (44) 

v,y  p 

where  [20] 

I  2 

Tij(v.u)  =  2v/Yi  |  fdR  n°iv*(R)  Waj (R )  n°ju(R)|  A(Ei v"Ejy,  Tjy)  (45) 

with 

Me.r)  =  n-1  (f/2)/(c2+f2/4)  (U6) 
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Here  1%^  is  the  half  width  of  level  p  in  WJj  due  to  phonon  transitions.  Because 
T(iv,ju)  describes  transitions  between  discrete  states  of  the  unperturbed  hamiltonian, 
care  must  be  exercized  to  include  the  width  of  the  initial  state.  This  is  done  for¬ 
mally  by  replacing  the  energy  conserving  6-function  by  the  Lorentzian  (46). 

Vie  will  now  report  on  numerical  results  on  thermal  field  evaporation  of  tungsten 
ions  from  a  tungsten  tip.  The  latter  we  model  by  a  finite  size  cluster  of  tungsten 
atoms,  which,  to  mimick  an  isolated  atom  on  a  (111)  surface,  is  chosen  as  three 
atoms  in  a  plane  with  one  atom  below  the  midpoint  of  the  triangle,  and  a  fifth  one 
in  the  symmetrical  position  above  the  plane.  Our  calculations  are  based  on  the 
semi-empirical  ASED-MO  method,  described  in  earlier  papers  [12,  18];  ionization  ener¬ 
gies  and  wavef unctions  for  W  are  those  of  Ref -[12]  except  that  we  have  raised,  in  an 
ad  hoc  manner,  the  ionization  energies  for  the  adatom  by  1.5eV  so  that  we  can 
describe  by  one  set  of  parameters  both  the  neutral  and  ionic  state  of  this  atom.  In 
this  approach  we  unfortunately  cannot  incorporate  the  electric  field  in  a  selfconsis- 
tent  manner.  However,  because  the  local  variation  of  the  electric  field  is  impor¬ 
tant,  we  take  the  electric  field  from  selfconsislent  calculations  [14,15]  for  a 
plane  jellium  metal  using  density  functional  theory  assigning  a  Wigner-Seitz  radius 
rs=2.07  to  tungsten.  We  then  impose  this  field  onto  the  cluster  assuming  that  the 
jellium  edge  is  half  a  lattice  spacing  above  the  plane  of  atoms  in  our  cluster. 

To  determine  the  transition  probabilities  (45)  we  must  know  the  nuclear  wave  func¬ 
tions,  n°ov  and  h°iui  in  the  diabatic  potentials,  W00  and  Wu.  To  simplify  the  numerics 
we  have  fitted  a  Morse  potential  to  W00,  i.e. 

W00(z)  =  A0(F)  [exp[-2Y(z-z0)  ]  -  2exp[-  Y(z-z0)  ]]  (47) 

adjusting  its  parameters,  A0(F),  Y(F)  and  z0(F),  as  a  function  of  field  strength. 
Likewise,  we  set 

Wu(z)  =  Wc  -eF(z-zc)  )48) 

for  the  diabatic  curve  of  the  ion.  For  both  potentials  the  wave  functions  can  be 
given  analytically. 

Field  evaporation  being  an  activated  process,  it  is  instructive  to  parametrize  the 
ion  yield  rate  constant  according  to  the  Polanyi-Wigner  equation  (1).  In  Table  I  we 
present  the  relevant  data.  We  note  first  that  the  activation  energies  for  field  eva¬ 
poration,  Q(F),  are  in  good  agreement  with  experiment.  To  put  this  into  perspective, 
we  want  to  point  out  that  this  agreement  depends  on  our  choice  for  rs,  as  discussed 
in  [1],e.g.  if  we  take  rs=1.5  the  calculated  evaporation  field  strength  increases  bya 
about  25%  from  its  value  at  rs=2.07. 

Turning  next  to  the  prefactor  v(F)  we  first  note  its  dependence  on  field  strength. 
In  the  early  days  of  field  ion  spectroscopy  it  had  been  assumed  that  the  prefactor 
is  independent  of  field  strength  [21].  However,  Kellogg  [7]  has  found  a  substantial 
field  dependence  in  field  evaporation  of  tungsten.  For  thermal  field  desorption  of 
helium  our  theory  [  1 8]  explained  such  a  field  1-  ;  i'-'.oe  as  the  result  of  the 
changes  in  the  shape  of  the  surface  potential  in  which  helium  is  adsorbed.  For  field 
evaporation  this  effect  is  also  present  albeit  of  lesser  significance.  To  obtain  a 
better  understanding,  we  recall  that  in  thermal  desorption  in  the  absence  of  a  field, 
the  desorption  rate  constant  at  low  coverage  can  be  written  for  desorption  from  a 
mobile,  nonlocalized  adsorbate  as 

rd  =  Svzexp(-Ed/kBT)  (49) 
when  hvz<<kgT,  and  as 


1 1 


rd  =  S  kBT/h  exp(-Ed/kBT)  (50) 

when  hvz>>kBT.  From  a  localized  adsorbate  we  get 

rd  =  S  2/ mas/kBT  vxvxvzrxp(-Ed/kgT)  (51) 


where  vx  etc.  are  the  vibrational  frequencies  of  the  adsorbed  particle  at  the  bottom 
of  the  surface  potential  well  in  the  x,  y,  and  z  directions.  Note  that  the  desorp¬ 
tion  rate  also  contains,  due  to  detailed  balance,  the  sticking  coefficient  S  as  a 
measure  of  the  efficiency  of  energy  transfer  between  the  solid  and  the  adatom. 

For  field  desorption  and  field  evaporation  the  adatom  must  be  thermally  excited  up 
to  the  energy  of  the  potential  barrier;  this  is  a  process  akin  to  thermal  desorption. 
Subsequently,  it  must  get  ionized.  Thus  the  prefactor  consists  of  two  factors, 
namely  an  ionization  probability  a(F)  and  an  attempt  frequency  v(F).  Their  field  dep¬ 
endences  are  opposing  each  other  in  that  a(F)  increases  from  zero  in  zero  field  to 
one  at  "high  "  fields,  whereas  v(F)  decreases.  We  can  view  this  process  as  a  parti¬ 
cle  localized  in  the  excited  level,  ei=Q(F),  at  the  energy  of  the  potential  barrier 
attempting  with  a  frequency  vi=ei/h  of  that  level  to  ionize  with  a  probability  a(F). 
Because  Q(F)  decreases  with  F,  so  does  and  thus  vj.  In  addition,  the  hump  of  the 
activation  barrier  and  thus  the  region  for  ionization  moves  towards  the  metal  for  in¬ 
creasing  field  resulting  in  an  increasing  ionization  probability  a(F).  In  Table  1  we 
have  estimated  the  critical  energy  level  in  W00,  asumed  to  be  a  Morse  potential,  and 
find  that  e.g.  for  F= ^ . 7 V/$  the  adatom  is  in  the  37th  excited  state  when  it  attempts 
to  ionize.  This  level  has  a  frequency  v37=2.86x1 O'-s-1  as  opposed  to  the  ground  state 
frequency  v„=4.5x1 012s_1.  This  compares  rather  well  with  the  prefactor  of 
2.76x1  O^s-1  obtained  from  the  Polanyi-Wigner  parametrization  of  the  ion  yield.  This 
argument  should,  however,  not  be  taken  too  literally  as  other  factors  contribute  to 
the  prefactor  as  one  already  knows  from  the  simpler  situation  of  thermal  desorption, 
cf.  equations  (49-51).  In  particular  this  estimate  does  not  allows  us  tc  extract 
a(F)  because  it  would  obviously  be  larger  than  one,  contrary  to  its  definition.  How¬ 
ever,  we  can  guess  that  for  fields  less  than  4.5 V/A  the  field  dependence  of  a(F) 
dominates,  and  most  likely  a(F)=1  for  larger  fields,  where  the  decrease  in  v(F) 
becomes  dominant. 

Although  our  theory  produces  the  right  trend  in  the  prefactor,  namely  decreasing 
with  increasing  field  for  fields  larger  than  4.5V/JI,  as  observed  by  Kellogg,  there 
are  discrepancies  in  the  absolute  values  in  that  the  experimental  data  are  substan  - 
tially  lower  except  at  F=4.7V/ft.  However,  Kellogg  warns  that  his  prefactors  have  an 
uncertainty  of  at  least  one  order  of  magnitude.  On  the  other  hand,  our  theoretical 
values  may  also  be  out  by  an  order  of  magnitude  due  to  several  of  our  approxima¬ 
tions,  in  particular  neglecting  lateral  variations  in  the  electric  field  which  might 
effect  the  localization  of  the  adatom.  Also  recall  that  the  master  equation  that 
underlies  (44)  as  derived  in  refs.  [I6,l8],is  restricted  in  its  applicability  to  the 
markovian  limit.  With  prefactors  of  the  same  order  (10l“s-1)  as  thermal  phonon 
assisted  transitions  in  W00,  one  should  account  for  nonmarkovian  effects  by,  e.g., 
using  the  Tsukada-Gortel  equations  [20], which  will  result  in  lower  prefactors,  but 
most  likely  not  more  than  one  order  of  magnitude. 

We  would  like  to  comment  further  on  the  difference  in  prefactors  for  field  evapo¬ 
ration  of  tungsten  as  opposed  to  thermal  field  desorption  of  helium.  For  the  latter 
case  we  found  [18]  prefactors  increasing  more  or  less  exponentially  from  a  low  of 
106s-1  at  F=4V/&  to  a  high  of  1012s_l  at  6V/S.  This  has  been  interpreted  as  due  to  an 
increase  in  the  ionization  probability  a(F)  due  to  rapidly  increasing  overlap  of  wave 
functions  on  the  helium  and  in  the  metal.  Afterall,  in  the  same  field  strength 
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regime  the  activation  barrier  for  helium  increases  substantially  whereas  for  metal 
field  evaporation  the  activation  barrier  decreases  monotonically ,  i.e.  for  helium  the 
field  enhances  adsorption,  at  least  up  to  6V/A.  Eventually,  this  enhancement  ceases, 
most  likely  around  7V/A  due  to  the  same  effects  that  cause  a  monotonic  weakening  of 
the  surface  bond  for  metals  on  metals,  namely  a  draining  of  bonding  orbitals  into 
the  metal. 

We  have  so  far  presented  results  on  thermal  field  evaporation  of  singly  charged 
tungsten  ions.  Experiment,  on  the  other  hand,  only  detects  W3+  and  W“+  ions.  The 
fact  that  our  calculated  prefactors  are  rather  large,  i.e.  of  the  order  of  10“*s-1, 
we  take  as  evidence  to  suggest  that  the  first  ionization  stage  from  W  to  W+  must  be 
the  slowest,  i.e.  rate  determining  step.  Because  the  tunneling  rate  into  the  metal 
decreases  rapidly  with  increasing  distance,  the  higher  ionization  states  cannot  be 
produced  by  tunneling  into  the  metal  as  their  abundance  would  in  that  case  be  decre¬ 
ased  over  singly  charged  ions  roughly  by  the  ratio  of  the  ionization  rates.  Thus 
post  ionization  [24]  happens  at  least  several  angstroms  away  from  the  surface,  but 
still  in  the  high  field  region,  by  field  ionization  with  the  excess  electrons  tunnel¬ 
ing  into  vacuum  states  rather  than  into  empty  metal  states.  Simple  minded  estimates 
confirm  this  scenario  although  detailed  calculations,  e.g.  for  tungsten,  are  not  ava¬ 
ilable  at  this  stage. 

6.  Conclusions 

In  this  paper  we  have  applied  our  recent  theory  of  the  kinetics  of  field  ionization, 
field  desorption  and  field  evaporation  to  study  field  evaporation  of  tungsten.  We 
have  developped  a  perturbative  method  to  calculate  diabatic  states  from  adiabatic 
ones.  We  find  good  agreement  of  the  field  dependence  of  the  activation  barrier  with 
experimental  data.  We  also  find  that  the  prefactor  in  the  ion  yield  increases  in¬ 
itially  (up  to  4.5V/A  for  tungsten)  due  to  an  rapid  increase  in  the  ionization  proba¬ 
bility  and  then  decreases  with  increasing  field  strength  due  to  changes  in  the  sur¬ 
face  potential.  There  is  some  discrepancy  in  the  absolute  value  which  must  be 
resolved  by  better  experiments  and  by  a  better  theory.  For  the  latter,  we  are  cur¬ 
rently  adopting  a  cluster  programme  based  on  the  spin  density  functional  theory  to 
more  properly  account  for  charged  species  by  a  better  treatment  of  Coulomb  effects. 
We  will  then  also  present  results  for  other  metals  to  hopefully  stimulate  further 
experiments . 
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Table  1 

Field  dependence  of  the  activation  barrier,  Q,  and  the  prefactor,  av,  experimental 
data  from  Ref. [7].  Also,  the  ground  state  frequency  v0  and  the  excited  frequency 
for  which  hvi=Q;  for  details  see  text. 


Field 

(V/A 

Q(ev) 

exp.  calc. 

k0(sec-1) 
exp.  calc. 

v0(sec-1) 

ni 

v^=(2n^+1 ) v 
(sec-1) 

3.0 

- 

1.85 

- 

1  .OIxlO1** 

2.92x1 012 

99 

4.58x10“* 

3.50 

- 

1.72 

- 

1.83x10“* 

4.93x1 012 

55 

5.47x1 011* 

n.o 

- 

1.42 

- 

2.33X1011* 

4.55x1 012 

49 

4.54x10“* 

4.5 

- 

1 .22 

- 

2.93x10“* 

4.62c 1 012 

41 

3.78x10“* 

4.70 

0.90 

1.17 

3x1 016 

2.76x1 011* 

4.52x1 012 

37 

2.86x10“* 

4.93 

0.60 

0.85 

1  xl  013 

2.60x10“* 

4.00x1 012 

30 

2.1 1x10“* 

5.10 

0.52 

0.58 

8x1 012 

2.21  xl  011* 

3.51x1 012 

24 

1 .58x10“* 

5.30 

0.35 

0.41 

7x1 0U 

2.25x1  O11* 

3.1  7x1 012 

18 

1  .lOxlO11* 

5.47 

0.31 

0.31 

7x10n 

2.06x10“* 

2.83x1 O12 

15 

8.30x1 013 

5.72 

0.20 

0.20 

3x1 011 

1.66x10“* 

2.68x1 O12 

10 

5.37x1013 

5.92 

0.12 

0.12 

4x1 011 

5.46x1013 

2.64x1 O12 

7 

3.82x1 013 

